ON THE PARAMETER DOMAIN OF WISHART DISTRIBUTIONS 
AND THEIR INFINITE DIVISIBILITY 
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Abstract. A complete characterization of Wishart distributions on the cones of positive 
scmi-dcfinitc matrices is provided in terms of a description of their maximal parameter 
domain. This result is new in that also degenerate scale parameters are included. For 
such cases, the standard constraints on the range of the shape parameter may be relaxed. 
Furthermore, the infinitely divisible Wishart distributions are revealed as suitable trans- 
formations and embeddings of one dimensional gamma distributions. This note completes 
the findings of Levy (1937) concerning infinite divisibility and Gindikin (1975) regarding 
the existence issue. 



1. Introduction 

The Wishart distribution T{p;a) on the cone S£ of symmetric positive semi-definite 
d x d matrices is defined (whenever it exists) by its Laplace transform 

C( 1 (p;a))(u) = (det(I + au)y p 1 UE-S+, (1.1) 

were p > denotes its shape parameter and a G S£ is the scale parameter. In the non- 
degenerate case where a is invertible, and for a discrete set of shape parameters, these 
distributions have been introduced in 1928 by Wishart [7] as sums of "squares" of centered 
multivariate normal distributions (quite similarly to the the chi-squared distributions in 
dimension one). 

In 1937, Levy jl] showed that a) on is not infinitely divisible for invertible a, 
which means that for some sequence of shape parameters pk I 0, 7(pfc;c) cannot exist. 
Gindikin [2] , Shanbhag pj and Peddada & Richards subsequently showed: 

Theorem G. For non- degenerate a the following are equivalent: 

(i) The right side of (II. ip is the Laplace transform of a probability measure. 

(ii) p belongs to the Gindikin ensemble 

A d = U, j = l,2,...,d-2|u 



Key words and phrases. Wishart distribution, Gindikin ensemble, matrix variate distributions, param- 
eter domain. 

1 Contrary to [5] we exclude the point mass at zero, i.e. the Gindikin ensemble does not contain 0. 
That's why we have chosen p > 0. 
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Aim of the present note is to extend this characterization by also allowing for degenerate 
a G Sj, and, to determine all Wishart distributions which are infinitely divisible. 

It is worth noting that - in view of the right side of (11. ip and Levy's continuity theorem 
- the maximal parameter domain 

@d := {(p, a) G K + x S£ \ T(p; a) is a probability measure} 

is closed, and that only for invertible scale parameters the issues of existence and infinite 
divisibility have been settled in the literature (where one easily follows from the other). 
Hence the following two statements answer very natural question on the Wishart family. 

Theorem 1.1. Let d G N, p > and a G S£ with rank(er) = r. The following are 
equivalent: 

(i) The right side of (II. ip is the Laplace transform of a non-trivial probability measure 
T(p;a) on Sj. 

(ii) p G A r . 

Theorem 1.2. Let d G N, p > and a G Sj[. The following are equivalent: 

(i) T(p; a) is infinitely divisible. 

(ii) rank(a) = 1. 

Remark 1.3. • It should be pointed out that Theorem 1 1 . 1 1 does not follow from the 

classification of positive Riesz distributions, when rank(cr) < d, see, e.g., [TJ Theo- 
rem VII. 3. 2]. In fact let us consider the special case, where a = diag(l, 0, . . . , 0) G 
Sj, and denote by 7id-^ r the projection onto the r-th subminors of Sj~, that is 

Kd-tr '■ S d -)■ S r , 7T d ^ r (a) = (dij)l<ij< r . (1.2) 

Then the right side of (II .ip takes the form 

(det(7r^ r (l + u)))-v = (1 + u u )- p , (1-3) 

which the Laplace transform of an infinitely divisible probability measure on St 
due to the theorems I1.1HTT21 However, the "corresponding" Riesz distribution □ 
(which is also infinitely divisible) has Laplace transform 

( d e t( ^( ( / + u) - 1))) P = (- i i±^_) P , 

which obviously differs from eq. (I1.3p . Nevertheless, both functions are charac- 
teristic functions of probability measures on the sub-cones of positive matrices of 
rank 1. 

• It is possible to generalize the Theorems of this paper to symmetric cones. In that 
more general setting the set A must be replaced by the so-called Wallach set, see, 
e.g., [TJ Theorem VII. 3.1]. We avoid this setting to make the presentation short 
and accessible to a larger group of readers. 

2 This is a member of the natural exponential family generated by the standard Riesz distribution with 
Laplace transform (det(7r ( i_ ) . r (w~ 1 ))) p 
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Notation: Id is the d x d unit matrix, and if there arises no confusion, we simply write 
I. det and tr are determinant and trace operator, and rank denotes the matrix rank. 
For matrices A, B of r x r and t x t dimension, diag(A, B) denotes the corresponding 
block-diagonal (r + t) x (r + t) matrix. 

Proof of Theorem \l.l[ Let U be an orthogonal matrix such that UaU T = diag(.D, 0, . . . , 0) 
where D = diag(<7i, . . . , a r ). We introduce the linear automorphism gu as the map 

gu : S+ S+, 9u (C) := £/^ T . 



Proof of (i) =>■ (ii) : We denote by T* the push- forward measure of T(p; a) under gt/, which 
means that for Borel sets A G Sj we have r»(A) = r(p; (x)(^ 1 (A)) = T(p; a)(U T AU). By 
[5], eq. (1.3.4)] the functional determinant of gu equals 1, hence for all s G —Sj[ we have 

£(T,)(«):=/ e te «>r*(dO= / e tr ^r(p; a)(d([/ T ^)) 



= det(/ + c7f/ T sf/)- p = det (J + f/af/ T s)- p 

= det(/ d + diag(£>, 0)s)~ p = det(J r + Ds)~ p . (1.4) 

By Theorem G, we have that there exists a probability measure T{p; D) on Sj[ with Laplace 
transform 

det(/ r + Du)~ p , u G -S+. 

We may therefore conclude that r* equals r(p; D)*' nd ^ r , the pullback of r(p; D) under 
Tid^r (see the defining equation (11.21) ). In fact, for all s G —S£ we obtain 

£(r(p; £>)*•*<-■)(*) = / e tr( ^ (s)€) r(p ;J D)((iO, 

M 

which equals eq. (II. 4p . But since diag(o"i, . . . , o~ r ) is of full rank r, we have by Theorem G 



that p G A n and we have shown that (ii) holds. 

Proof of (ii) =>- (i) : We reverse the preceding arguments: By Theorem G, there exists a 
probability measure T{p, D) on for p G A r . Let ir r ^ d be the embedding 

TT r ^ d : S+ -)• S+, 7r r ^ d (A) := diag(A 0). 

Then the Laplace transform of the push-forward T(p, D)1 r ^ d of T(p, -D) under n r ^ d equals 

£(r(p, L>)f^) = det(/ d + diag(L>, 0, . . . , 0)7^(5))-* 

Hence by eq. (II. 4p the push- forward of T(p, D)1 r ^ d under g v -\ = g^ 1 = gu T (hence given 
by % 1 (0 = U T t;U) equals T(p,a), which proves (i) □ 



Proof of Theorem \1.2l (i) (ii) : By infinite divisibility and the right side of eq. (II. ip we 
see that T(p/n; a) is a probability measure for each n > 1. But by definition, C [1/2, oo) 
if and only if d > 1. Hence by Theorem ll.il we have that rank(a) = 1. 



4 



EBERHARD MAYERHOFER 



The converse implication (ii) => (i) may be proved in a similar way: By Theorem 11.11 and 
in view of the fact that Ti = (0, oo), we have that T(p/n; a) is a probability measure for all 
n > 1. Therefore by the right side of eq. (11. ip we see that T(p; a) is infinitely divisible. □ 
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